Phase transitions in the Shastry-Sutherland lattice 
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Two recently developed theoretical approaches are applied to the Shastry-Sutherland lattice, vary- 
ing the ratio J' / J between the couplings on the square lattice and on the oblique bonds. A self- 
consistent perturbation, starting from either Ising or plaquette bond singlets, supports the existence 
of an intermediate phase between the dimer phase and the Ising phase. This existence is confirmed 
by the results of a renormalized excitonic method. This method, which satisfactorily reproduces 
the singlet triplet gap in the dimer phase, confirms the existence of a gapped phase in the interval 
0.66 < J'/J < 0.86 



I. INTRODUCTION 

Strongly correlated electron systems, and in partic- 
ular those which can be treated as quantum magnets, 
are the subject of a continuous and intense attention 
from experimentalists and theoreticians. They frequently 
exhibit complex phase diagrams, and non-trivial low en- 
ergy physics. Among the spin lattices those present- 
ing spin frustrations both present unconventional phases 
and quantum phase transitions and represent a challenge 
for theoreticians. The quantum Monte Carlo methods, 
which may be considered as the most rehable treatment 
for non-frustrated 2-D or 3-D lattices, cannot be ap- 
plied in these cases. Among the 2-D frustrated lattices 
one may quote the checkerboard, the Kagome and the 
triangular lattices, which have been the subject of nu- 
merous theoretical works. The SrCu2{BO^)2 lattice^ 
is a famous two-dimensional anti-ferromagnetic system 
presenting a spin gap and free from long range order. 
The Copper atoms are of (fi character and can be seen 
as S' = 1/2 spins. This lattice may be considered as 
a realization of the Shastry-Sutherland model^ which 
can be schematized as a square lattice, with J' anti- 
ferromagnetic coupling between nearest neighbors, and 
diagonal anti-ferromagnetic interactions J in one plaque- 
tte over two, as pictured in Fig. ^ This system is sup- 
posed to obey the corresponding Heisenberg Hamiltonian 



H = J' ^ ^ SiSj + J ^ ^ SiSj 



(1) 



where the couples [i,j] concern the bonds of the square 
lattice and the connected pairs of next nearest 

neighbor atoms. The real material has been the sub- 
ject of intense experimental studies, showing the exis- 
tence of a spin gap^ an almost localized nature of the 
singlet triplet excitation;^ and the existence of magneti- 
zation plateaux»^ The Shastry-Sutherland Hamiltonian 
has been widely studied by theoreticians (for review see 
S. Miyahara and K. Ueda-), varying the J'/J ratio. For 
small values of J'/ J the ground state is a product of 
singlets strictly localized on the oblique bonds. The 
real material would correspond to a ratio J'/J — 0.635, 
close to the critical value where this phase disappears. 



Oppositely when J is small, the perturbed 2-D square 
lattice is gapless with long range order. Early stud- 
ies based on exact diagonaHzations^ Ising expansion^ 
or dimer expansions and real space renormalization 
group with effective interaction (RSRG-EI}li predict a 
simple transition between the Ising phase and the dimer 
phase for J'/J = 0.696. Other works have suggested 
the existence of an intermediate phase. From Schwinger 
boson mean field theorj*iS this phase would be an heli- 
cal ordered state, ranging between J'/J = 0.6 and 0.9 
while a plaquette expansioni^*!^ predicts a plaquette sin- 
glet based phase for 0.677 < J'/J < 0.86, a result sup- 
ported by other exact diagonalization results Weihong 
et alpifi suggested that the intermediate phase between 
J'/J = 0.69 and 0.83 might be columnar rather than 
plaquette singlet. 

The present paper employs two recently developed 
methods to study the phase diagram of this lattice. We 
first will consider the cohesive energy of the lattice, us- 
ing a self consistent pertubation (SCP)ii which may be 
seen as a modified coupled cluster expansion li^iiSiSSiSi 
The method starts from a reference function which can 
be the Ising spin distribution or a columnar product of 
bond singlets. It happens that the cohesive energy calcu- 
lated from the latter reference function is the lower one 
for 0.66 < J'/J < 0.86, the lowest energy being obtained 
from the Ising reference for J'/J > 0.86. The calculated 
cohesive energies obtained by the RSRG-EI with (3 x 3) 
sites square blocks or by the contractor renormalization 
(CORE)Si with columnar blocks of (2 x 6) sites agree 
with the values obtained from SCP and support the idea 
of an intermediate phase. 

A second section concentrates on the calculation of 
the gap, using a renormalized excitonic method (REM), 
based on the scale-change ideas inspiring the renormal- 
ization group methods. The methods proceeds through 
the definition of blocks with even number of sites. The 
design of the relevant blocks is different for the dimer 
phase and for the other phases. The calculated gap for 
the dimer phase is in good agreement with that of pre- 
vious theoretical studjii2*i^ and compatible with experi- 
mental evidencesii Approaching the problem from colum- 
nar blocks one obtains a finite gap for J'/J < 0.88. This 
result confirms the existence of an intermediate phase. 
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II. CALCULATION OF THE COHESIVE 
ENERGY 

If one defines A = J'/ (J' + J) and take J' + J as the 
energy unit, the product of singlets on the oblique dimers 
(J interactions) is an eigenfunction whatever A, and its 
cohesive energy is J' + J = 1. 



A. Self-consistent perturbation 

We have used the recently proposed SCP methodii 
to evaluate the cohesive energy of other phases. The 
method can be seen as a modified coupled cluster 
methodpiSiiSiSSi21 It starts from a reference function $0; 
supposed to be a relevant zero-order function for the con- 
sidered phase. This function is highly locaHzed. In prac- 
tice $0 will be 

- the Ising spin distribution on the square 2-D lattice, 

- or a product of bond singlets in a columnar arrange- 
ment. On each of these bonds one may also define 
a local triplet state. 

The action of H on $0 generates a first generation of 
local excited function $i (($i|i/|$o) 7^ 0). In the so- 
called intermediate normalization convention 



* = *o + I] ^^^^ + Yl 



(2) 



where the vectors <i>Q do not interact with $0 
{{^a\H\<^Q) = 0), the knowledge of the coefficients of 
the first generation $i's is sufficient to fix the ground 
state energy 



$o|i^l$^ 



(3) 



Starting from Ising <i>o, the vectors <i>i are obtained by 
a spin exchange on any bond of the 2-D square lattice 
(Fig. nj. Starting from the product $q of singlets on 
the X-directed bonds in a columnar arrangement (Fig. 
I2II, the vectors <i>i are excited singlet states produced as 
products of two triplets on interacting bonds. They are 
of four different types in a Shastry-Sutherland lattice, as 
pictured in Fig. [2l 

The determination of the coefficients d is governed by 
the corresponding eigenequation, adopting the compact 
notation i/y = ($i|i?|<i>j). 

An estimation of the coefficients of the second generation 
functions $0, (such that Hia)is only necessary for those 
which are obtained from by operations in the strict 
neighborhood of the bonds involved in the process 
creating $i from $0 (*&* = T^^o)- The spin exchanges 
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FIG. 1: Genealogy of the wave function starting from Ising 
"I'd. *&o = reference function, $i = first generation function, 
"I'd, $/3 types of non factorisable second generation func- 
tions. 



or excitations on remote bonds cancel, according to 
the linked cluster theorem, since in that case 



(5) 



The number of second generation functions to be con- 
sidered is therefore very limited. Their coefficient is 
estimated from the C^'s according to perturbative ar- 
guments, with a systematic consideration of exclusion 
principle violating corrections of the energy denomina- 
tors, which introduce infinite summations of diagrams 
and speed the convergence. Details of the method are 
given elsewhercfli The equations for the precise problem 
are given in the Appendix, the results appear in Fig. El 
One sees that 

- the energy calculated from Ising $0 is the lowest 
one for large values of J' /J. This energy curve 
cuts the energy of the locaHzed dimer phase for 
J' /J = 0.69. This value compares quite well with 
other estimates (SiSiiSiii and in particular with that 
of a RSRG-Elii with (3 x 3) sites square blocks 
{[J'/J)c = 0.696) 

- however the energy calculated from the columnar 
arrangement of bond singlets on the 2-D square lat- 
tice happens to be lower than the preceding one for 
J'/ J < 0.859. The corresponding energy is the low- 
est one for 0.661 < J'/ J < 0.859. Hence these cal- 
culations support the suggestion of the existence of 
three phases, with an intermediate phase between 
the dimer phase and the Ising phase. 
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FIG. 2: Genealogy of the wave function starting from colum- 
nar product of bond singlets empty ellipses. Dark ellipses 
represent local triplets. The figure pictures the 4 types of 
first generation singlet states. The non factorisable second 
generation functions are too numerous to be pictured. 



However the fact that one finds two distinct values of 
the energy from two distinct reference functions within 
an approximate algorithm does not prove the existence 
of two phases. It is possible that at convergence the two 
wave operators, fl from $0 and fl' from $q, lead to the 
same wave function 



As an argument in favor of two distinct phases we may 
mention that starting from an other function $q, ener- 
getically degenerate with % |$[,') = m\H\%)), 
product of bond singlets on parallel bonds which do not 
belong to the same plaquette (see Fig. 0J, the calcu- 
lated energy always remains above that obtained from 
Ising $0, as seen in Fig. El It is not easy, at least from 
our SCP approach, to determine wether the intermedi- 
ate phase is based on a columnar arrangement of bond 
singlets or is a plaquette phase. The second one keeps 
isotropic properties in the X and Y directions, while the 
former one breaks this symmetry. It is interesting to 
remark that our reference wave function $q is strongly 
anisotropic, with a probability zero to find the parallel 
spins in the X-directed bonds of a plaquette and a prob- 
ability 1/2 to find parallel spins in the F-directed bonds. 
After consideration of the possible excitations, the proba- 
bility to find parallel spins along the two AT-directed bond 
of an empty plaquette becomes 87% of the probability to 
find parallel spins on the F-directed bonds of the same 
plaquette. The isotropy is almost restaured by the self- 
consistent evaluation of the excitation amplitudes. The 
result, obtained from a highly broken-symmetry reference 
function, pleads in favor of a plaquette phase. 



FIG. 4: Alternative reference function $0'. 
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FIG. 3: Cohesive energy (with the unit {J + J')) as a function 
of the A parameter. The horizontal straight line concerns the 

dimer phase, SCP expansion from columnar bond singlets 

"3?0i ■ ■ ■ SCP expansion from Ising 'E'o, — SCP expansion from 
shifted bond singlets $0'. 



B. Renormalization group methods 

Renormalization group techniques have also been em- 
ployed. A preceding work^^ had considered square 9 sites 
i.e., (3 X 3) sites blocks, the ground state of which is 
an 5*2 — ±1/2 doublet, i.e., a quasi spin. The lattice 
of blocks is isomorphic to the Shastry-Sutherland lat- 
tice, with two types of effective interactions J'^^^ and 
J*^^) between the blocks, the amplitude of which can 
be obtained from the exact spectrum of the dimers of 
blocks, according to the effective Hamiltonian theory.— 
The ratio J''-^^ / is larger than J'/ J for J' /J > 0.696. 
For J'/ J = 0.696, J'(i)/j(i) = J'/ J, i.e., this value is 
a fixed point, for which the cohesive energy calculated 
by iterating the RSRG-EI procedure practically coincide 
with that of the dimer phase. The behaviour of the so- 
calculated cohesive energy of the Ising phase is plotted 
in Fig. El The intersection of the cohesive energy curves 
of the dimer phase and of the RSRG-EI estimate of the 
Ising phase occurs for J' / J = 0.672iil The method does 
not produce an intermediate phase, when it is applied to 
these isotropic square blocks. 
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TABLE I: Cohesive energy (with the unit (J + J')) for inter- 
mediate values of the J' j 3 ratio. 



J' /J 


SCP(Ising) 


SCP(Plaquette) 


CORE*^ 


RSRG-El'3 


0.65016 


-0.97804 


-0.99257 


-0.99803 


-0.98501 


0.66667 


-0.98987 


-1.00355 


-1.00787 


-0.99609 


0.69461 


-1.00989 


-1.02191 


-1.02488 


-1.01490 


0.72413 


-1.03026 


-1.04034 


-1.04251 


-1.03362 


0.81818 


-1.09291 


-1.09604 


-1.09783 


-1.09510 


0.85185 


-1.11419 


-1.11475 


-1.11720 


-1.11790 


1.0 


-1.20061 


-1.19028 


-1.19478 


-1.20920 



(a) from (2 x 6) sites columnar blocks 

(b) from (3 x 3) sites square blocks (cf. ref. 11) 



If one uses 2n sites blocks, with a non degenerate sin- 
glet ground state of energy E\ , the ground state of 
the lattice is studied from the product of the ground 
state of the blocks = Ha V'a according to the sim- 
plest version of the CORE methodi^ The knowledge of 
the ground state exact energy E\g of the AB dimers en- 
ables one to define an effective interaction vab between 
the blocks^ as 



VAB 



El 



El 



and the energy per block is 

Ea = EI + \Y.''ab- 



(7) 



(8) 



In the research of a tentative columnar phase, we have 
considered 12 (2 x 6) sites columnar blocks built from 
three aligned plaquettes. Interestingly enough the cal- 
culated cohesive energy is the lowest one (cf Fig. in 
the interval 0.656 < J' / J < 0.845 (RSRG-EI) or 0.901 
(SCP). This result is in excellent agreement with our SCP 
result, and supports the existence of an intermediate pla- 
quette phase. Table 1 reports the calculated values of the 
cohesive energy in the problematic domain of the J' / J 
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FIG. 6: Evolution of the cohesive energy (with the unit (J -I- 
J')) as a fonction of A = J' /{J + J')- The horizontal line 
concerns the dimer phase, - - - SCP from columnar bond 
singlets <E>o, ■ ■ ■ SCP from Ising <E>o, RSRG-EI from (3x3) 
sites square blocksJA □ CORE with (2 x 6) sites columnar 
blocks. The lower part is a zoom in the region of crossing 
between the intermediate and Ising phases. 



ratio one sees the good agreement between the SCP from 
Ising and RSRG-EI on one hand and between the SCP 
from columnar arrangement of bond singlets and CORE 
with (2 X 6) sites columnar blocks on the other hand. 




(A) (B) 

FIG. 5: (A) definition of the periodizable 12 sites blocks for 
the calculation of the excitation gap in the dimer phase, (B) 
definition of (2 x 2n) sites columnar blocks for the calculation 
of the excitation gap in the plaquette phase, and application 
of the CORE method. 



III. CALCULATION OF THE SINGLET 
TRIPLET GAP 

The study of the singlet triplet gap may bring addi- 
tional informations. It has been performed according to 
an other renormaHzation group technique, the REM(2i 
which starts from blocks A with non degenerate ground 
state singlet and an excited triplet state V'^) ^^"^ build 
the excited state from linear combinations of locally ex- 
cited states, ^\ = V'aIIbV'b- The knowledge of the 
exact spectrum of the dimers of blocks makes possible 
the calculation of the effective interaction between 



an excited triplet and neighbor singlets and of an effec- 
tive excitation-hopping integral Hab which couples V'aV'b 
with V'aV'b- 

We first have applied the method to the dimer phase. 
It is known that in this phase the triplet excitation on 
a given dimer bond can only propagate to a neighbor 
one through a 6th-order process involving a vortex of 4 
dimers. A complete evaluation of the interaction leads to 

a 



{SabTcd\H-ff^'^\TabS,d) = 



2. J' 



/6 



J5 




d 
-• 



This interaction remains small when J' increases, and the 
energy gap slowly deviates from J. A convenient design 
of 12 sites (6 dimers bonds) blocks involving two vor- 
tex have been pictured in Fig. EJA). The gap calculated 
from the REM are reported in the left part of Fig. to- 
gether with those of a dimer expansioni^ The excitation 
energy is tangent to the line A/ J = 1 at the origin, as ex- 
pected, but falls down rapidly near the phase transition. 
It vanishes for J'/J = 0.7006. For the J'/ J = 0.635, i.e., 
the value proposed for the real material, our calculation 
gives A — 0.32 J, which satisfactorily compares with the 
experimental value (30-35 K)i if one accepts the usually 
proposed value of J (85 K). For the columnar plaque- 
tte phase, blocks built from aligned plaquettes have been 
considered (Fig. ElB)), together with an extrapolation 
on the block sizei^i The results appear in Fig. ^together 
with those of a plaquette expansion One sees that 
the system is gapped for J'/J < 0.883, and gapless be- 
yond this value. The gap in this phase goes through a 
maximum for J'/ J ~ 0.73, a result also obtained by Koga 
and Kawakamifi^ii^ Our results confirm the existence of 
an intermediate gapped phase. 




IV. CONCLUSION 



The present work, employing essentially different 
methods based on either coupled cluster type expansions 
or renormalization group techniques, present consistent 
results confirming the existence of an intermediate phase 
in the Shastry-Sutherland lattice, and its plaquette or 
columnar nature. The results regarding the location of 
the phase transitions are gathered in Table 2. While the 
phase transition between the dimer and the Ising phases 
would occur at J'/J between 0.67 and 0.70, the interme- 
diate phase is the lowest one from J'/J = 0.656 (CORE 
12 sites) - 0.661 (SCP) to J'/J = 0.859 (SCP) - 0.901 
(CORE/SCP). Obtained from new methods, these re- 
sults agree with those of refs. 13-16 and reduce the re- 
maining uncertainties concerning the existence, domain 
and nature of the intermediate phase in the Shastry- 
Sutherland lattice. 



TABLE II: Critical values of the J/ J' ratio, based on cohesive 
energy (Coh) or gap calculations (Gap). 





Criterion 


Methods 


J'/J 


Dimer/Ising 


Coh 
Coh 


dimer'^/SCP'^ 
dimer^/RSRG-El'^ 


0.690 
0.672 




Gap 


REM'^ 


0.701 


Dimer/Interm. 


Coh 
Coh 


dimer^/SCP^ 
dimer*/CORE^ 


0.661 
0.656 




Coh 


SCP^/SCP'' 


0.859 


Interm. /Ising 


Coh 
Coh 


CORE^SCP^ 
SCP*'/RSRG-EI'= 


0.901 
0.826 




Coh 
Gap 


coreVrsrg-ei<= 
remS 


0.845 
0.883 



(a) exact energy 

(b) from Ising $0 

(c) from (3 x 3) sites square blocks (cf. ref. 11) 

(d) gap vanishing of the dimer phase, 12-site blocks 
(Fig. IIA)) 

(e) from columnar bond singlets 'I'o 

(f) from (2 X 6) sites columnar blocks 

(g) gap vanishing of the Intermediate phase, 
extrapolations from columnar blocks. 



FIG. 7: Dependence of the singlet-triplet gap on the J'/J 
ratio. — REM from 12-site blocks (Fig. ElA)) for the dimer 

phase, results of a dimer expansion;^ o extrapolated REM 

from columnar plaquette blocks (Fig. ElB)), □ results of a 
plaquette expansion^i— 
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APPENDIX A 



SCP equations for the Shastry-Sutherland lattice 

(1) Equation from Ising function. There is a unique type of spin exchange, and a single coefScient C. 



[6 J' - 2 J - 7 J'C] C + J' + JC + J' 
"l2J'-4J- 14 J'C 



+4J' 



10 J' - 2 J - 13 J'C 
0, 



+ 10 J' 



12J'-4J- 14 J'C \ _ 
8 J' - 4 J - 12 J'C ) ~ 

12 J' - 4 J - 14 J'C 



10 J' - 4 J - 13 J'C 



1 



+ 2 J' 



/ 12J'-4J-14J'C \ _ ^ 



8J'-2J- 12 J'C / 

12 J' - 4 J - 14 J'C 
12 J' - 2 J - 14 J'C ~ 



8 J' - 2 J - 8 J'C 



(Al) 



Ecoh = -2J'(1 - C). 



(A2) 



(2) Equations from bond singlets functions. There are four types of double excitations, as pictured in Fig. 2. Using 
the notations 



the four coupled polynomial equations are 

[2 J' - J - {hiCi + 2/12C2 + 4/i3C3 + 2/14C4)] Ci + hi + ""'"3 



2hiCi 



2/12 (/12C1 + hiC2) + 4/i3(2/i3Ci + 2/iiC3) + 2/14 (/i4Ci + /11C4) 



4J' - (2/iiCi + 2/12C2 + 4/13C3 + 2h4Ci) 



0, 



(A3) 



r , / XT 2 rft,4C3 + /13C4] 2/l2C 

[2 J' - (2/iiCi + /12C2 + 4/13C3 + 2/14C4)] C2 + /i2 - ^ ^^ ^ + ^ ■ 



V3 



2/li(/tiC2 + /l2Ci) + 2/t3(2/l3C2 + 2/t2C3; 

4J' - (2/iiCi + 2/12C2 + 4/13C3 + 2/14C4) 



0, 



(A4) 



[3 J' - (2/.1C1 + 2/.2C2 + 3h,C, + 2h,C,)] C3 + /.3 - + ^^[^' + ^' + ^4] 



\/3 ■ 3 

2fei(2/tiC3 + 2/l3Ci) + h2{2h2C3 + 2/13C2) + 4fe^C3 + fe4(2/l4C3 + 2/t3C4) 

4J' - (2/iiCi + 2/12C2 + 4/13C3 + 2/14C4) 



= 0, 



(A5) 



r , . M 2 r/l2C3 + /l3C2] 2/14 C| 

4j' - J - (2/iiCi + 2/12C2 + 4/13C3 + /14C4)] C4 + /i4 - ^ V + -^r^ 

v3 3 



2/11 (/11C4 + hiCi) + 2h^{2h^Ci + 2/14 C3) 
4J' - (2/iiCi + 2/12C2 + 4/13C3 + 2/14C4) 



= 0, 



{A6) 



[ - (7 J' + J) + /iiCi + /12C2 + 2/13C3 + /14C4] 



(A7) 
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